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I.  INTRODUCTION 


The  accurate  prediction  of  projectile  aerodynamics  is  of  significant 
importance  in  the  early  design  stage  of  a  projectile.  An  area  of  considerable 
interest  in  transonic  projectile  aerodynamics  is  the  determination  of  the 
critical  aerodynamic  behavior.  The  critical  aerodynamic  behavior  of  projec¬ 
tiles  occurs  in  the  transonic  speed  regime,  0.9  <  M  <  i.l,  where  the  aerody¬ 
namic  coefficients  have  been  found  to  increase  by  as  much  as  100%.  This  rapid 
change  in  the  aerodynamic  coefficients  can  be  attributed  in  part  to  the  com¬ 
plex  flow  structure  existing  on  projectiles  at  transonic  speeds.  For  projec¬ 
tiles  flying  at  angle  of  attack,  the  shock  pattern  is  asymmetric  and  thus, 
contributes  to  the  critical  aerodynamic  behavior. 

In  recent  years  a  considerable  research  effort  has  been  focused  on  the 
development  of  modern  predictive  capabilities  for  determining  projectile 
aerodynamics.  Numerical  capabilities  have  recently  been  developed  using  the 
Navi er-Stokes  computat i onal  technique  to  compute  the  flow  over  slender  bodies 
of  revolution  at  transonic  speeds.  This  technique  was  used  to  predict  pro¬ 
jectile  aerodynamics  for  both  axi symmetri c 1  and  3-D  flow2  situations.  The  3-D 
flow  predictions2  were  rather  poor  which  could  be  due  to  lack  of  sufficient 
grid  resolution  and,  more  importantly,  lack  of  correct  base  flow  modeling. 
These  calculations  did  not  include  the  base  region  flow.  The  upstream  effect 
the  base  region  flow  may  have  on  the  afterbody  flow  field  and  the  associated 
shock  structure  was,  thus,  neglected.  Deiwert2  has  also  successfully  computed 
the  three  dimensional  flow  field  over  a  boattailed  afterbody.  Emphasis  in  his 
work  was  on  the  boattail  flow  field  and  not  on  the  base  region  flow.  Recent¬ 
ly,  an  axisymmetric  thin-layer  Navier-Stokes  base  flow  code  has  been  developed 
to  compute  the  entire  projectile  flow  field  including  the  base  region  using  a 
unique  flow  field  segmentation  procedure.4  This  code  has  been  used  to  compute 
base  flow  with  and  without  mass  injection4  and  including  jet  effects. 5  The 
objective  of  this  initial  effort  is  to  extend  the  capability  of  this  code  to 
three-dimensions  in  order  to  include  the  upstream  effect  of  the  base  flow  and 
to  calculate  the  pitch  plane  projectile  aerodynamics  at  transonic  speeds.  The 
asymmetric  shock  pattern  existing  on  the  projectile  and  possible  base  region 
influence  are  felt  to  be  the  major  contributors  to  the  critical  aerodynamic 
behavior  at  transonic  speeds. 

This  report  describes  the  initial  capability  which  is  used  to  predict  the 
flow  field  over  a  6  caliber  secant-ogi ve-cyl inder-boattai  1  projectile  at  M  = 
.96  and  a  =  4°.  Computed  results,  obtained  on  a  Cray  X-MP  computer  are  pre¬ 
sented  and  compared  to  experimental  data.  A  brief  description  of  the  govern¬ 
ing  equations  and  the  flow  field  segmentation  procedure  is  included. 


II.  GOVERNING  EQUATIONS 

The  Navier-Stokes  equations  solved  herein  make  use  of  the  thin-layer 
approximation.6  That  is,  the  viscous  terms  involving  velocity  gradients  in 
both  the  longitudinal  and  circumferential  directions  are  neglected.  The 
viscous  terms  are  retained,  however,  for  velocity  gradients  in  a  direction 
nearly  normal  to  the  surface  where  large  flow  field  gradients  exist.  This 
formulation  retains  the  momentum  equations  in  all  three  coordinate  directions. 
The  retention  of  the  three  momentum  equations  allows  for  the  computation  of 
separated  flow  and  thus  differs  s i gn 1 f i cant ly  from  boundary-1 ayer  assumptions. 
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The  equations  solved  here  are  written  in  a  generalized  coordinate  system. 
This  allows  a  wide  variety  of  body  shapes  to  be  computed  using  the  same  basic 
numerical  technique.  The  three-dimensional,  transformed,  thin-layer  Navier 
Stokes  equations  written  in  nondimensional ,  strong  conservation  law  form  are7 


'6  q  +  3,E  +  3  F  +  3  G  =  Re-^3,S 


n  t 


where  generalized  coordinates 


£,  =  £(x,y,z,t)  -  longitudinal  coordinate 
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The  velocities  in  the  £,  n,  c  coordinates  are 


u  =  h  +  Cxu  +  V  *  *zw 

V  =  nt  +  nxu  +  nyv  +  nzw  (2) 

w  =  ;t  +  cxu  +  SV  +  ;zw 

which  represent  the  contravari ant  velocity  components. 

The  Cartesian  velocity  components  (u,  v,  w)  are  retained  as  the  dependent 

variables  and  are  nondimensi onal i zed  with  respect  to  a^  (the  free  stream  speed 

of  sound).  The  local  pressure  is  determined  using  the  relation 


p  =  (y  -  l)[e  -  .5p(u2  +  v2  +  w2)]  (3) 

where  y  is  the  ratio  of  specific  heats,  density  (p),  is  referenced  to  p^ 
and  total  energy  (e)  to  P^a^2. 


III.  NUMERICAL  METHOD 


1.  COMPUTATIONAL  ALGORITHM 

The  numerical  scheme  used  for  the  solution  of  Equation  (1)  is  a  fully, 
implicit,  approximately  factored,  finite  difference  algorithm  in  delta  form  as 
analyzed  by  Beam  and  Warming.8  This  scheme  can  be  first-  or  second-order 
accurate  in  time  and  second-  or  fourth-order  accurate  in  the  three  spatial 
directions.  The  solution  of  the  three-dimensional  equations  is  implemented  by 
an  approximate  factori zati on  which  allows  the  system  of  equations  to  be  solved 
in  three  one-dimensional  steps.  This  procedure  has  been  utilized  in  previous 
applications2'9  with  high  degree  of  success.  Additional  details  of  the  numer¬ 
ical  method,  computati onal  algorithm  and  boundary  conditions  can  be  found  in 
Reference  9. 
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2.  FLOW  FIELD  SEGMENTATION 


A  unique  flow  field  segmentation  procedure  has  been  used  for  three- 
dimensional  flow  over  a  projectile  including  the  base  flow.  This  is  similar 
to  the  one  previously  used  in  the  axisymmetric  flow  situation4  where  two  2-D 
adjoining  grids  were  used.  The  segmentation  process  allowed  the  existing 
finite  difference  algorithm  to  be  applied  to  two  3-D  adjoining  grids  by  making 
modifications  in  the  internal  structure  of  the  block  tridiagonal  matrix  and 
boundary  conditions.  This  greatly  simplified  the  development  of  the  code  for 
computations  of  flow  over  a  projectile  at  angle  of  attack  including  the  base 
regi  on. 


Figure  1  is  a  schematic  illustration  of  the  flow  field  segmentation  used 
to  compute  the  entire  projectile  flow  field  that  includes  the  base  flow.  It 
shows  h  w  the  two  grids  A8C0  and  AEFG  in  the  physical  plane  are  transformed 
into  a  single  segmented  grid  in  the  computati ona 1  domain.  An  important  advan¬ 
tage  of  this  segmentation  procedure  lies  in  the  preservation  of  the  sharp 
corner  at  the  base,  which  allows  easy  blending  of  the  computational  meshes 
between  the  regions  ABCD  and  AEFG.  No  approximation  of  the  actual  sharp 
corner  at  the  base  is  made.  Thus,  realistic  representation  of  the  base  is 
inherent  in  the  current  procedure.  The  cross-hatched  region  represents  the 
projectile.  The  line  BC  is  the  projectile  base  and  the  region  ABCD  is  the 
base  or  wake  region.  The  line  AB  is  a  computational  cut  through  the  physical 
base  region  that  acts  as  a  repetitive  boundary  in  the  computational  domai  n. 


With  the  computational  domain  so  segmented,  implicit  integration  is  car¬ 
ried  out  in  c,  t,  and  n  directions  by  solving  the  block  tridiagonal  matrices 
first  in  C  and  then  in  ;  and  n.  The  block  tridiagonal  matrix  in  the 
4  direction  has  elements  at  J  =  JB,  J8+1  which  are  treated  as  internal  boun¬ 
daries  in  the  computati onal  domain  (J  =  JB  represents  the  projectile  base  and 
J  =  J 8+1  is  the  nose  axis).  The  block  tridiagonal  matrix  in  the  E,  direction 
takes  the  following  form  (after  setting  Lj  =  0  to  simplify  the  illustration) 


1  1  U  p 
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•AJB+1  ! 


JMAA-2 


■'JB 
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—  ~ 
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• 
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• 
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• 
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0 
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• 

• 
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• 
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• 

AqJMAX-l 

RHSJMAX-1 
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Here  A 1 s  denote  the  quantity  -^-7  A ,  I  is  a  5  *  5  identity  matrix  and  RHS  is 
t re  right  hund  side  of  the  finite-difference  form  of  Equation  (1).  Note  the 


appearance  of  the  uncoupled  block  tri di  agona  1  S. 
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The  rows  at  JB  and  JB+1  are 


•HP 


particularly  simple  because  boundary  conditions  are  updated  explicitly  at  the 
end  of  inversion*.  One  simply  fills  the  block  tridiagonal  matrix  ignoring  the 
base  JB  and  the  nose  axis  JB+1.  Elements  in  these  rows  are  then  overloaded  as 
shown  above.  The  flow  field  segmentation  does  not  affect  the  block  tridiagon¬ 
al  matrix  in  the  ;  and  n  directions.  This  is  readily  accomplished  by  initial¬ 
ly  ignoring  these  boundaries  and  then  overwriting  the  tridiagonal  elements. 


3.  BOUNDARY  CONDITIONS  AND  TURBULENCE  MODEL 

The  no  slip  boundary  condition  for  viscous  flow  is  enforced  by  setting  the 
contravari ant  velocities  to  zero,  i.e. , 


U  =  V  =  W  =  0 


on  the  projectile  surface.  The  flow  field  in  the  immediate  near-wake  region 
has  been  considered  to  be  weakly  viscous  and  inviscid  boundary  condition  is 
used  at  the  base.  The  viscous  terms  near  the  projectile  surface  and  across 
the  shear  layer  are  considered  to  be  the  most  do'.l( r:a;i'  and  have  been  retained 
within  the  thin-layer  approximation  used  here.  This  approximation  neglects 
the  viscous  terms  normal  to  the  base  that  are  believed  to  be  less  significant. 

Recently,  solutions  have  been  obtained  for  axisymmetric  base  flow  computa¬ 
tions  where  viscous  terms  in  the  longitudinal  direction  or  normal  to  the  base 
were  included.  These  unpublished  results  obtained  with  full  Navi er-Stckes 
equations  and  viscous  boundary  condition  at  the  base  are  not  found  to  be  sig¬ 
nificantly  different  from  the  earlier  solutions  obtained  with  the  thin-layer 
approximation  and  inviscid  base  boundary  condition.  This  justifies  our  use  of 
inviscid  boundary  condition  at  the  base.  Additionally,  the  grid  points  used 
for  the  three  dimensional  calculation  results  in  a  rather  coarse  grid  in  the 
near  wake.  The  pressure,  p  is  calculated  using  the  three  transformed  momentum 
equations  and  density,  p  is  extrapolated  at  the  inner  boundary  (body  surface). 
At  the  base,  these  two  quantities  are  extrapolated. 

At  the  outer  boundary,  constant  free  stream  values  are  used  for  all  vari¬ 
ables.  The  flow  variables  above  and  below  the  cut  (line  AB  in  Figure  1)  are 
simply  averaged  to  determine  the  boundary  conditions  along  the  cut.  All 
variables  at  the  downstream,  upstream  axis  and  on  the  center  line  of  the  wake 
region  are  obtained  through  extrapolation.  For  nonspinning  projectiles  at 
angle  of  attack,  symmetry  exists  about  the  axis  of  the  projectile  and  the 
computation  is  performed  only  over  the  half  plane.  Thus,  symmetry  boundary 
conditions  are  imposed  at  both  ends  in  the  ci rcumferenti al  direction.  Because 
of  this  boundary  condition,  only  the  static  aerodynamic  coefficients  for  pitch 
plane  aerodynamics  can  be  determined.  For  spinning  projectiles  at  angle  of 
attack,  the  full  three  dimensional  flow  field  must  be  computed. 

For  the  computation  of  turbulent  flows,  a  turbulence  model  must  be  sup¬ 
plied.  In  the  present  calculations,  a  two-layer  algebraic  eddy  viscosity 
model  due  to  Baldwin  and  Lomax6  is  used.  In  their  two-layer  model,  the  inner 
region  follows  the  Prandtl-Van  Driest  formulation.  The  outer  formulation  can 
be  used  in  attached  and  separated  boundary  layers  and  with  minor  modifications 
in  wakes  as  well.  In  both  the  inner  and  outer  formulations,  the  distribution 
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of  vorticity  is  used  to  determine  the  length  scales,  thereby  avoiding  the 
necessity  of  finding  the  outer  edge  of  the  boundary  layer  (or  wake).  For  base 
or  wake  flow,  a  few  minor  modifications  are  made.  These  are  discussed  in 
Reference  4  and  are  not  included  here.  The  algebraic  eddy  viscosity  model  is 
simple  to  use;  however,  it  is  not  strictly  valid  for  such  wake  flows.  Two 
equation  turbulence  models  are  better  suited  for  these  flows.  These  models 
have  been  used  to  compute  axi symmetric  base  flow10  and  have  failed  to  provide 
better  predictions  than  those  obtained  with  simple  algebraic  eddy  viscosity 
model.  imp  re  ..cnt  in  the  advanced  two-equation  turbulence  modeling  is 
needed.  The  accuracy  of  the  numerical  computations  should  improve  with 
improvement  in  the  turbulence  modeling. 


IV.  MODEL  GEOMETRY  AND  EXPERIMENT 

Une  means  of  establishing  the  computational  accuracy  of  a  numerical  scheme 
is  through  comparisons  with  available  experimental  data.  The  model  used  for 
the  experiment  and  computati ona 1  study  presented  here  is  an  idealization  of  a 
realistic  artillery  projectile  geometry.  The  experimental  model  shown  in 
Figure  2  consists  of  a  three-cal i  ber  (l  caliber  =  maximum  body  diameter), 
sharp,  secant-ogive  nose,  a  two-caliber  cylindrical  midsection,  and  a  one- 
caliber  7-deg  conical  afterbody  or  boattail.  A  similar  model  was  used  for  the 
computati onal  studies  with  the  only  difference  being  a  5%  rounding  of  the  nose 
tip.  The  nose  tip  rounding  was  done  for  computational  efficiency  and  is  con¬ 
sidered  to  have  little  impact  on  the  final  integrated  forces. 

T^e  experimental  pressure  data11  12  used  for  comparison  in  this  report 
were  obtained  in  the  NASA  Langley  8-ft  Pressure  Tunnel  using  a  sti  ng-inounted 
model.  The  test  conditions  of  1  atm  supply  pressure  and  320  K  supply  tempera¬ 
ture  resulted  in  a  Reynolds  number  of  4,5  *  106  based  on  model  length.  The 
pressure  data  were  obtained  with  the  model  nonspinning. 


V.  RESULTS 

As  indicated  earlier,  this  solution  technique  involves  solving  the  time- 
dependent,  thin-layer,  Na vi er-Stokes  equations.  The  procedure  is  started  by 
assuming  uniform  free-stream  conditions  for  all  grid  points  in  the  computa¬ 
tional  domain.  A  slow  start  of  the  boundary  conditions  is  implemented  and  the 
calculation  marches  in  time  until  a  steady-state  solution  is  obtained.  A 
criterion  for  convergence  is  for  the  solution  residual  to  decrease  by  three 
orders  of  magnitude.  Also,  the  surface  pressure  distribution  is  checked  for 
time  invariance.  The  implicit  technique  used  here  allows  for  large  computa¬ 
tional  time  steps  to  be  taken,  which  helps  to  reduce  the  total  computation 
time.  All  the  computed  results  to  be  presented  are  for  M  =  .96  and  a  =  4°. 

One  of  the  first  steps  in  performing  a  computation  is  the  generation  of  a 
computational  grid.  Tnese  points  are  determined  prior  to  tne  computations  and 
are  not  changed  with  time.  An  example  of  the  computational  grid  used  in  this 
study  is  shown  in  Figure  3.  This  figure  shows  m  expanded  view  of  the  grid 
for  the  boattailed  afterbody  and  the  base  region.  The  surface  of  the  boattail 
and  the  base  can  be  seen  along  with  the  grids  in  the  windward  and  leeward 
planes.  These  grids  (excluding  the  base  region)  were  generated  using  a  modi¬ 
fied  version  of  an  elliptic  grid  generation  code  described  in  Reference  13. 

6 


The  grid  in  the  base  region  was  then  added  following  the  flowfield  segmenta¬ 
tion  procedure  to  obtain  the  full  grid.  The  outer  flowfield  boundary  was 
placed  approximately  18  body  diameters  from  the  body  surface.  At  this  dis¬ 
tance,  the  flowfield  should  be  uniform  and  the  imposed  free-stream  boundary 
conditions  are  valid.  The  clustering  of  grid  points  near  the  body  surface  is 
required  to  resolve  the  viscous  boundary  layer  near  the  body  surface.  Grid 
clustering  has  also  been  used  in  the  longitudinal  direction  near  the  boattail 
and  the  base  corners  where  large  gradients  in  the  flow  variables  are  expected. 
The  largest  grid  size  that  could  be  accommodated  in  the  Cray  X-MP/22  computer 
at  NASA  Ames  Research  Center  was  70  X  36  X  21  (longitudinal,  normal  and  cir¬ 
cumferential).  This  grid  resolution  is  inadequate  for  numerical  simulation  of 
the  entire  projectile  including  the  base  region.  In  order  to  obtain  meaning¬ 
ful  results,  computations  were  carried  out  in  two  pieces.  The  first  one 
corresponds  to  the  solution  over  the  projectile  with  emphasis  on  the  forebody 
solution  X/D  <  4.5  and  the  other,  X/D  >  4.5,  containing  the  boattail  corner 
and  the  base  region.  For  each  of  these  cases,  a  70  X  36  X  21  grid  was  used. 
Computations  were  first  run  for  the  projectile  where  the  base  flow  was  modeled 
as  an  extended  sting.2  The  converged  solution  at  X/D  =  4.6  was  then  used  as 
an  upstream  boundary  condition  for  computation  of  the  boattailed  afterbody  and 
base  region  flow  field  using  the  grid  shown  in  Figure  3.  The  code,  however, 
has  been  developed  so  that  complete  numerical  simulation  of  a  projectile 
including  the  base  region  can  he  made.  Such  simulations  require  larger  com¬ 
puter  memory  and  are  currently  in  progress  using  the  BRL  CRAY  X-MP/48 
computer. 

Particle  paths  in  the  base  region  for  both  the  windward  and  the  leeward 
planes  are  shown  in  Figure  4.  This  figure  shows  the  general  feature  of  the 
flow  in  the  base  region.  The  reci rcul  atory  flow  in  the  base  is  evident.  In 
addition,  this  figure  shows  the  expected  asymmetry  in  the  reci rculatory  flow 
pattern  in  the  wind-side  and  the  lee-side  planes.  The  size  of  the  separation 
bubble  in  the  lee-side  is  much  larger  than  that  on  the  wind-side.  Figure  5 
shows  the  grid  and  the  Mach  number  contours  in  the  same  perspective  view  which 
includes  the  windward  and  leeward  planes  as  \  e 1 1  as  the  circumferential  plane 
at  the  end  of  the  afterbody,  i.e.  ,  at  the  base.  As  seen  in  Figure  5a,  the 
grid  stretching  at  the  base  region  in  the  normal  direction  is  used  on  both 
sides  of  the  base  comer.  The  corresponding  Mach  contours  in  Figure  5b  show 
the  expansion  wave  front  in  the  base  region  which  seems  to  grow  as  the  lee- 
side  is  approached. 

Figures  6  and  7  are  two-dimensional  views  which  show  the  Mach  contours  and 
pressure  contours,  respectively,  over  the  boattail  and  in  the  base  region  for 
the  windward  and  leeward  planes.  Figure  6  shows  mixed  regions  of  subsonic 
flow  and  supersonic  flow.  The  recirculation  region  in  the  near  wake  is  sub¬ 
sonic.  It  also  indicates  the  asymmetry  in  the  flow  field  in  the  base  region 
and  the  presence  of  strong  free  shear  layers.  Pressure  contours  are  shown  in 
Figure  7  which  show  the  expansion  at  the  boattail  corner  and  indicate  the  pre¬ 
sence  of  a  shock  wave  on  the  boattail.  The  grid  in  the  longitudinal  direction 
on  the  boattail  is  rather  coarse  and,  thus,  a  sharp  shock  wave  is  not  seen. 
The  expected  asymmetry  in  the  near  wake  flow  field  is  clearly  seen  in  this 
figure.  The  position  of  the  shock  wave  is  better  shown  in  Figure  8.  The 
snoefc.  location  is  based  cn  looking  at  the  Mach  number  component  in  the  direc¬ 
tion  of  the  local  pressure  gradient.  The  shock  structure  is  asymmetric.  The 
snock  wave  on  the  boattail  in  the  '-'ind-side  is  closer  to  the  base  than  the  one 
in  the  lee-side. 
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The  surface  pressure  coefficient  is  shown  as  a  function  of  the  longitudi¬ 
nal  position  for  both  the  wind-side  and  the  lee-side  in  Fiqure  9.  The  compu¬ 
tational  results  are  compared  with  the  experimental  data1"*  shown  in  squares 
and  circles.  The  agreement  between  the  computed  surface  pressure  coefficient 
and  experimental  data  on  the  wind-side  is  good.  The  expansions  and  recompres¬ 
sions  near  the  ogive-cylinder  and  cyl i nder-boattai 1  junctions  are  captured 
adequately  by  the  computation.  Comparison  of  the  surface  pressure  on  the  lee- 
side  shows  discrepancy  on  the  ogive  as  well  as  the  boattail.  The  agreement 
with  experimental  data  is  better  near  the  ogi ve-cyl i nder  junction.  The  com¬ 
puted  surface  pressure  on  the  lee-side  is  lower  than  that  on  the  wind-side; 
whereas  the  reverse  is  true  on  the  boattail.  A  similar  trend  is  observed 
experimental ly.  Comparisons  between  computation  and  experiment  for  the 
ci rcumferential  surface  pressure  distribution  ere  shown  in  Figures  10-13  for 
longitudinal  positions  X/D  s  4.5,  5.19,  5.56  and  5.78,  respectively.  Figure 
10  shows  the  comparison  at  X/D  =  4.5.  This  is  on  the  cylindrical  portion  up¬ 
stream  of  the  boattail.  The  agreement  between  the  computation  and  experiment 
is  excellent.  This  solution  was  used  as  an  upstream  boundary  condition  for 
the  boattailed  afterbody  and  base  region  flow  computation.  Comparisons  at  the 
next  three  stations  are  on  the  boattail.  Figures  11  and  12  show  the  compari¬ 
sons  at  X/D  =  5.19  and  5.56,  respecti vely.  As  seen  in  these  figures,  the  com¬ 
putation  predicts  the  correct  trend  of  the  data.  The  agreement  is  good  near 
the  wind-side  (<f>  =  0)  and  gets  worse  as  the  lee-side  (if  =  100°)  is  approached 
(overpredicting  at  X/D  =  5.19  and  underpredicting  at  X/D  =  5.56).  The  next 
station,  X/D  =  5,78  is  on  the  boattail  just  ahead  of  the  base  corner  (Figure 
13)  and  here  the  influence  of  the  base  region  flow  should  be  the  largest. 
Discrepancy  between  the  computed  result  and  the  experiment  exists  at  all 
points  in  the  ci rcumferential  direction  at  this  station;  however,  the  trend  of 
data  is  again  correctly  captured  in  the  computed  result.  The  discrepancy  is 
partly  due  to  the  fact  that  the  experimental  data  were  obtained  using  a  sting 
mounted  model  whereas  the  computation  does  not  include  the  sting  and  also 
partly  due  to  lack  of  adequate  grid  resolution  on  the  boattail.  The  effect 
that  the  asymmetrical  ly  located  shock  structure  existing  on  the  boattail  can 

have  on  the  pitching  moment  is  not  clearly  understood  and  is  a  subject  of 

further  investigation. 

The  base  pressure  for  this  S0CBT  projectile  has  been  obtained  and  plotted 
as  a  function  of  roll  angle  in  Figure  14.  The  computed  result  (solid  line)  is 
compared  to  experimental  date1^  shown  in  circles.  This  base  pressure  is 
taken  at  44%  of  the  base  radius  from  the  center  line.  As  seen  in  this  figure, 
t^e  base  pressure  coefficient  seem  to  decrease  slightly  in  going  from  the 
wind-side  to  the  lee-side.  This  trend  is  also  observed  in  the  experiment. 
The  computed  result  is  underpredicted;  however,  it  should  be  noted  that  these 
base  pressures  were  measured  with  a  sting  mounted  in  the  base  region.  The 

effect  of  the  sting  on  the  base  pressure  and  the  boattail  flow  field  has  not 
been  modeled  in  the  present  computation.  In  addition,  the  computational  grid 
’n  the  base  region  has  been  severely  stretched  in  the  streamwise  direction. 

The  final  result  desired  from  the  flowfield  calculations  is  the  determina¬ 
tion  of  the  aerodynamic  coefficients.  The  results  to  be  presented  here 

include  the  axial  force  and  normal  force  coefficients.  These  are  obtained  by 
integration  of  the  pressure  and  viscous  forces  acting  on  the  projectile. 
These  aerodynamic  coefficients  are  plotted  in  Figures  15  and  16  as  a  function 
of  longitudinal  position  and  illustrate  the  development  of  the  force  over  the 
length  of  the  projectile.  Fiqure  15  shows  the  axial  force  coefficient  as  it 
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begins  to  rise  over  the  first  half  of  the  nose  (P  >  PJ  and  then  drop  as  the 

ogive-cylinder  junction  is  approached  (P  <  Pj).  This  is  followed  by  a  sma1 

increase  over  the  cylinder  which  results  from  the  viscous  contribution  only. 
There  is  no  contribution  to  this  increase  due  to  the  pressure  forces  on  the 
cylinder.  The  axial  force  coefficient  rises  sharply  over  the  boattail  and 
then  levels  off.  Comparison  with  the  accumulated  axial  force  coefficient 
obtained  experimentally  shows  favorable  agreement.  The  experimental  result 
does  not  include  the  viscous  component;  and,  therfore,  is  expected  to  be  smal¬ 
ler  than  the  computati onal  result.  Figure  16  is  a  plot  of  the  normal  force 
coefficient.  The  rapid  increase  in  normal  force  that  occurs  on  the  ogive 
portion  of  the  projectile  is  shown.  The  cylinder  portion  should  produce  no 
significant  additional  normal  force;  however,  the  computation  indicates  a 
slight  increase  in  the  normal  force.  The  reversal  in  the  direction  of  the 
force  on  the  boattail  can  be  seen  clearly.  The  accumulated  normal  force  coef¬ 
ficient,  indicated  by  the  dark  dot,  has  been  cverpredi cted  by  the  computation. 
It  is  clear  that  the  boattail  has  a  dramatic  effect  on  the  aerodynamic  coeffi¬ 
cients.  Thus,  good  resolution  of  the  boattail  flowfield  is  essential. 


VI.  CONCLUDING  REMARKS 

This  report  described  the  development  and  the  application  of  an  unsteady, 
thin-layer.  Navi er- Stokes  base  flow  code  to  compute  the  3-D  base  flow  of 
projectiles  at  transonic  speeds. 

This  capability,  which  utilizes  an  implicit  finite  difference  scheme,  was 
used  to  predict  the  base  region  flowfield  of  a  projectile  at  Mm  =  .96  and  a  = 

4°.  The  solution  was  marched  in  time  for  approximately  2800  time  steps  to 
obtain  the  steady  state  result.  Computations  were  performed  on  a  Cray  X-MP/22 
computer  using  52920  grid  points  which  required  14  hours  computer  time. 
Computed  results  show  the  recirculation  region  and  the  asymmetric  wake  in  the 
base  region.  The  initial  results  indicate  the  presence  of  asymmetrically 
located  shock  structure  on  the  boattail  upstream  of  the  base  corner.  The  poor 
agreement  between  the  computed  surface  pressures  and  aerodynamic  coefficients 
with  experimental  data  indicate  the  need  for  more  grid  resolution.  These 
initial  results  are  encouraging  and  the  capability  has  now  been  established 
for  three  dimensional  numerical  simulation  of  a  projectile  including  the  base 
region.  The  need  for  additional  grid  resolution  exists  for  accurate  numerical 
predictions.  Accurate  prediction  of  all  these  coefficients  is  needed. 
Further  computations  are  then  needed  at  other  Mach  numbers  in  the  transonic 
reyirne  to  determine  the  critical  aerodynamic  behavior  of  projectiles. 


Figure  2.  Details  of  model  geometry 


Figure  4.  Particle  paths  in  windward  and  leeward  planes  in  the  base  region, 
=  .96,  u  =  4°. 
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Figure  14.  Circumferential  base  pressure  distribution. 
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LIST  OF  SYMBOLS 


maximum  speed  of  sound 

cross  sectional  area 

specific  heat  at  constant  pressure 

axial  force  coefficient  excluding  base  drag 

pitching  moment  coefficient 

normal  force  coefficient 

pressure  coefficient,  2(p-pJ/paou2 

body  diameter  (57,15mm) 

total  energy  per  unit  volume/p^a^ 

flux  vector  of  transformed  Navi er-Stokes  equations 

Jacobian  of  transformation 

Mach  number 

pressure/ Pooa£ 

Prandtl  number, u  c  /< 
vector  of  dependent  variables 
body  radius 

Reynolds  number,  p  a  0/p 
viscous  flux  vector 
physical  time 

Cartesian  velocity  comnone 'ts/a 

oo 

Contravari ant  velocity  component s/aw 
physical  Cartesian  coordinates 
angle  of  attack 
ratio  of  specific  heats 
coefficient  of  thermal  conductivity 
coefficient  of  viscosity 


29 


No.  of 
Copies 

12 


1 

1 


4 


1 


1 


1 


1 


DISTRIBUTION  LIST 
No.  of 

Organi zati on  Copies  Organi zati on 


Admi ni strator 

Defense  Technical  Info  Center 
ATTN:  DTIC-FDAC 
Cameron  Station,  Bldg  5 
Alexandria,  VA  22304-6145 

HQDA 

DAMA-ART-M 

Washington,  DC  20310 
Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

Commander 

Armament  RD&E  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-MSI 

SMCAR-LCA-F/Klein 

Loeb 

Hudgins 

Dover,  NJ  07801-5001 
Commander 

US  Army  Armament,  Munitions 
&  Chemical  Command 
ATTN:  SMCAR- IMP-L 
Rock  Island,  II  61299-7300 

rommander 

US  Army  Armament,  Research, 
Development  &  Engineering  Ctr 
ATTN:  SMCAR-TDC 
Dover,  NJ  07801 

Commander 

US  Army  Jefferson  Proving  Ground 
Materiel  Testing  Directorate 
ATTN:  Arthur  B.  Alphin,  MAJ ,  ARM 
Madison,  IN  47250-5100 

Commander 

CECOM  R&D  Technical  Library 

ATTN:  AMSEL-IM-L 

Fort  Monmouth,  NJ  07703-5000 


1  Director 

US  AMCCOM  ARDEC  CCAC 
Benet  Weapons  Laboratory 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  12189-4050 

1  Commander 

US  Army  Aviation  Systems  Cmd 

ATTN:  AMSAV-E 

4300  Goodfel low  Blvd. 

St.  Louis,  MO  63120-1798 

1  Director 

US  Army  Aviation  Research  & 
Technology  Activity 
Ames  Research  Center 
Moffett  Field,  CA  94035-1099 

3  Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-RX 

Redstone  Arsenal,  AL  35898-5249 

1  Director 

US  Army  Missile  &  Space 
Intelligence  Center 
ATTN:  AIAMS-YDL 
Redstone  Arsenal,  AL  35898-5500 

1  Commander 

US  Army  Tank  Automotive  Command 
ATTN:  AMST'A-TSL 
Warren,  MI  48397-5500 

1  Director 

US  Army  TRADOC  Analysis  Center 

ATTN:  ATOR-TSL 

White  Sands  Missile  Range, 

NM  88002-5502 

1  Commander 

US  Army  Research  Office 
P.  0.  Box  12211 
Research  Triangle  Park, 

NC  27709 


31 


No.  of 
Copies 


DISTRIBUTION  LIST 


No.  of 

Organi zation  Copi es  Organization 


1  Commander 

US  Naval  Air  Systems  Command 
ATTN:  AIR-604 
Washington,  DC  20360 

2  Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  Dr.  T.  Clare,  Code  DK20 
Or.  F.  Moore 

Dahlgren,  VA  22448-5000 

]  Commander 

US  Naval  Weapons  Center 
ATTN:  Code  3431,  Tech  Lib 
China  Lake,  CA  93565 

1  Commander 

US  Army  Development  & 
Employment  Agency 
ATTN:  MODE-ORO 
Ft.  Lewis,  WA  98433-5000 

1  Director 

NASA  Langley  Research  Center 
ATTN:  NS-185 ,  Teen  Lib 
Langley  Station 
Hampton,  VA  23365 

4  Director 

NASA  Ames  Research  Center 
ATTN:  MS-202- 1/Pul  1 i am 
MS-258-l/Steger 
Schi f f 
Jettmar 

Moffett  Field,  CA  94035 

2  Commandant 

US  Army  Infantry  School 
AT'N :  ATSH-CD-CS-OR 
Ft,  Benning,  GA  31905 

1  Commandant 
USAFAS 

ATTN:  ATSF-TSM-CN 

Ft.  Sill,  OK  73503-5600 


1  AFATL/DOIL  (Tech  Info  Ctr) 

Eglin  AFB,  FL  32542-5438 

2  Sandia  Laboratories 
ATTN:  Dr.  W.L.  Oberkampf 

Dr.  F.  Blottner 
Division  1636 

Sandia  National  Laboratories 
Albuquerque,  NM  87185 

1  AEDC 

Calspan  Field  Services 
ATTN:  MS  600  (Dr.  John  Benek) 
AAFS,  TN  37389 

1  Virginia  Polytechnic  Institute 
&  State  University 
ATTN:  Dr.  Clark  H.  Lewis 
Department  of  Aerospace  &  Ocean 
Engineering 
Blacksburg,  VA  24061 

1  University  of  California, 

Davi  s 

Department  of  Mechanical 
Engineering 

ATTN:  Prof.  H.A.  Dwyer 
Davis,  CA  95616 

1  Pennsylvania  State  University 
Department  of  Aerospace 
Engi neeri ng 

ATTN:  Dr.  G.  S.  Dulikravich 
University  Park,  PA  16802 

1  University  of  Florida 

Dept,  of  Engineering  Sciences 
College  of  Engineering 
ATTN:  Prof.  C.  C.  Hsu 
Gainesville,  FL  32611 

1U  CIA 

OIR/DB/Standard 
GE47  HQS 

Washington,  DC  20505 


1  AFWL/SUL 

Ki rtl and  AFB,  NM  87117 


32 


No.  of 
Copies 


DISTRIBUTION  LIST 


Organization 


1  University  of  Illinois 
at  Urbana  Champaign 
Department  of  Mechanical  $ 
Industrial  Engineering 
ATTN:  Prof.  W.  L.  Chow 
Urbana.  IL  61801 

1  University  of  Maryland 

Department  of  Aerospace  Engineering 
ATTN:  Dr.  J.D.  Anderson,  Jr. 
College  Park,  MD  20742 

1  University  of  Notre  Dame 

Department  of  Aeronautical 
&  Mechanical  Engineering 
ATTN:  Prof.  T.J.  Mueller 
Notre  Dame,  IN  46556 

1  University  of  Texas 

Department  of  Aerospace  Engineering 
&  Engineering  Mechanics 
ATTN:  Dr.  D.S.  Dolling 

Austin,  Texas  78712-1055 


K' 


Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP ,  H.  Cohen 

Cdr,  USATECOM 

ATTN:  AMSTE-SI-F 

Cdr,  CRDEC,  AMCCOM 
ATTN:  SMCCR-RSP-A 
SMCCR-M'J 
SMCCR-SPS-IL 


33 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Yo.  r  comments/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1 .  BRL  Report  Number _ Date  of  Report _ 

2.  Date  Report  Received _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report,  being  used? 
data,  proc  'ire,  source  of  ideas,  etc.) _ 


(Information  source,  design 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved 
etc?  If  so,  please  elaborate. _ _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


Name 


CURRENT 

ADDRESS 


Organization 


Address 


City,  State,  Zip 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD 

ADDRESS 


Organi r  at  ion 

Address 

City,  State,  Zip 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


